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Abstract 

In this talk we present the update (including the recent SNO results) of 
our analysis of the neutrino oscillation solutions of the solar and atmospheric 
neutrino problems in the framework of four-neutrino mixing where a sterile 
neutrino is added to the three standard ones and the mass spectra present two 
separated doublets. Such scenarios allow for simultaneous transitions of solar 
v e , as well as atmospheric Ua, into active and sterile neutrinos controlled by 
the additional mixing angles $23 an d $24 j and they contain as limiting cases 
the pure solar f e -active and v e -sterile neutrino oscillations, and the pure 
atmospheric — ► v s and — > v T oscillations, respectively. We evaluate the 
allowed active-sterile admixture in both solar and atmospheric oscillations 
from the combined analysis. Our results show that, although the Super- 
Kamiokande data disfavour both the pure Vn — > v s atmospheric channel and, 
in combination with SNO, the pure v e — * v s solar channel, the result from 
the combined analysis still favours close-to-pure active and sterile oscillations 
and disfavours oscillations into a near-maximal active-sterile admixture. 
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I. INTRODUCTION 



Super-Kamiokande high statistics data ||T]-0| indicate that the observed deficit in the 
/i-like atmospheric events is due to the neutrinos arriving at the detector at large zenith 
angles, strongly suggestive of the oscillation hypothesis. Similarly, their data on the 
zenith angle dependence and recoil energy spectrum of solar neutrinos in combination 



with the results from Homestake g], SAGE 0, and GALLEX+GNO |Ij|nj experiments 
and the recent SNO results ||12|| , have put on a firm observational basis the long-standing 
problem of solar neutrinos, strongly indicating the need for v e conversions. 

In addition to the solar and atmospheric neutrino results from underground experiments, 
there is also an indication for neutrino oscillations in the —>■ v e channel by the LSND 



experiment [I3|. All these results can be accommodated in a single oscillation framework 
only if there are at least three different scales of neutrino mass-squared differences. The 
simplest case in which this condition is satisfied requires the existence of a fourth light 
neutrino, which must be sterile (i.e. have interactions with standard model particles much 
weaker than the SM weak interaction) in order not to affect the invisible Z decay width, 
precisely measured at LEP [p^ -Tl8[|. 

There are various analyses in the literature that compare the oscillation channels into 
active or sterile neutrinos for both solar f2"2]-|2l)|j and atmospheric data [0,0 performed in 
the framework of two-neutrino oscillations (where oscillations occur into either an active or 
a sterile state). In particular the Super-Kamiokande collaboration has claimed that the os- 
cillation into sterile neutrinos is disfavoured for both atmospheric J7j and solar neutrinos || . 
The comparison of the recent SNO CC event rates with the Super-Kamiokande event rate 
for solar neutrinos imposes further constraints on this channel. However, when considered in 
the framework of four-neutrino mixing, oscillations into pure active or pure sterile states are 
only limiting cases of the most general possibility of oscillations into an admixture of active 



and sterile neutrinos ||18|| . Such a more general framework of four-neutrino mixing has also 
been analysed in recent studies, which have put constraints on the active-sterile admixture 
either in the context of solar |2f| or atmospheric |28H30| neutrino oscillations, respectively. 

In this talk we present the update of our combined analysis of the neutrino oscillation 
solutions of both the solar and the atmospheric neutrino problem, in the framework of 
four- neutrino mixing We include in our analysis the most recent solar neutrino rates 
of Homestake |], SAGE 0, GALLEX and GNO [0,0, as well as the recent 1258-day 
Super-Kamiokande data sample ||, including the recoil electron energy spectrum for both 



day and night periods and the recent results from the CC even rates at SNO |T2|]. As for 
atmospheric neutrinos we include in our analysis all the contained events from the latest 
79.5 kton-yr Super-Kamiokande data set |J, as well as the upward-going neutrino-induced 
muon fluxes from both Super-Kamiokande and the MACRO detector |J2], The constraints 



arising from the relevant reactor (mainly from Bugey and accelerator (CDHSW 
and CCFR ||35|| ) experiments are also imposed. 

The outline of the paper is the following. In Sec. |TJ we summarize the main expressions 
for the neutrino oscillation probabilities that we use in the analysis of solar and atmospheric 
neutrino data, which take into account matter effects both in the propagation in the Sun 
and in the Earth. Sections [TTI| and [IV] contain the results for the analysis of the four-neutrino 
oscillation parameters for solar and atmospheric data respectively. The results for the full 
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combined analysis are described in Sec. [V], with particular emphasis on the active-sterile 
oscillation admixture, which better allows for a simultaneous description of the solar and 
atmospheric neutrino data. Our results show that, against what may be naively expected, 
although the Super-Kamiokande data disfavour both the pure — > v s atmospheric channel 
and, in combination with SNO, the pure v e —>■ v s solar channel, the result from the combined 
analysis still favours any of these close-to-pure active and sterile oscillations and disfavours 
oscillations into a near-maximal active-sterile admixture. 



II. FOUR-NEUTRINO OSCILLATIONS 

There are six possible four-neutrino schemes that can accommodate the results from 
solar and atmospheric neutrino experiments as well as the LSND evidence. They can be 
divided in two classes: 3 + 1 and 2 + 2. In the 3 + 1 schemes there is a group of three 
neutrino masses separated from an isolated one by a gap of the order of 1 eV 2 , which is 
responsible for the short-baseline oscillations observed in the LSND experiment. In 2 + 2 
schemes there are two pairs of close masses separated by the LSND gap. The 3 + 1 schemes 



are disfavoured by experimental data with respect to the 2 + 2 schemes |fTE|,|PT| , but they are 
still marginally allowed ||36|| . Therefore we concentrate on the 2 + 2 schemes shown in Fig. [l|. 
For the phenomenology of neutrino oscillations these two schemes are equivalent up to the 
relabelling of the mass eigenstates (or equivalently of the mixing angles). Thus in what 
follows we will consider the scheme B, where the mass spectrum presents the hierarchy: 

Am 2 & = Am 2 ,! < Am 2 tm = Am 2 3 < Am 2 SND = Am 2 Al ~ Am 2 2 ~ Am 2 31 ~ Awf 2 (2.1) 

(we use the common notation Am^ = m| — m 2 ). In this four-neutrino mixing scheme, the 
flavour neutrino fields v a [a = e, s, fi, r) are related to the mass eigenstates fields v k by 

4 

v * = U °<k v k (a = e, s, n, r) , (2.2) 
fe=i 

where U is a 4x4 unitary mixing matrix. Neglecting possible CP phases, the matrix U can 
be written as a product of six rotations, U\ 2) Ui$, Uu, U23, U24 and U34 where 

(Uij) ab = 5 ab + (cij - 1) (S ia S ib + 6j a 8 jb ) + Sij (S ia 5 jb - 5 ja S ib ) , (2.3) 

where c^- = cos^y and = sin^. The order of the product of the matrices corresponds 
to a specific parametrization of the mixing matrix U. In order to study oscillations of the 
solar and atmospheric neutrinos, which include the matter effects in the Sun and/or in the 
Earth, it is convenient to use the following parametrization [T8| : 



U = U 24 U 23 U 1A U l3 U u U 12 . (2.4) 
This general form can be further simplified by taking into account the negative results from 



the reactor experiments, in particular the Bugey one [33], which implies that 



P e B e U9ev = l-2(\U el \ 2 + \U e2 \ 2 ) (l - \U el \ 2 - \U e2 \ 2 ) (sin 2 Am | SNpL ) > 0.99 (2.5) 
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in the range of Am^ relevant to the LSND experiment. This leads to an upper bound on 
the projection of the v e over the 3-4 states: 

l^esl 2 + |tfe4| 2 = 44 + 4 £ (2.6) 

so that the two angles $13 and $14 must be small and their contribution to solar and at- 
mospheric neutrino transitions is negligible. Therefore, in what follows we will set these 
two angles to zero, although it should be kept in mind that at least one of them must be 
non- vanishing in order to account for the LSND observation. 

In the approximation ^13 = $14 = the U matrix takes the form: 



U 



( Ci 2 S12 \ 

-S12C23C24 C12C23C24 S23C24C34 - S24S34 S23C24S34 + S 24C 3 4 

S12S23 -C12S23 C23C34 C23S34 

V S12C23S24 — C 12 c 23 s 24 — •S23S24C34 ~~ C 24-S34 — •S23S24S34 + C24C34 / 



(2.7) 



In the rest of this section we will discuss the consequences of the mixing structure of Eq. ( |2.7|) 
on the relevant transition probabilities for solar and atmospheric neutrino oscillations. The 



transition probabilities that take into account matter effects have been derived in Ref. |T8 
Some improvement concerning the calculation of the regeneration of solar f e 's in the Earth 
was presented in Ref. PB| . Here we summarize the discussion on the flavour composition of 
the relevant states for both solar and atmospheric neutrino oscillations. 



A. Solar Neutrinos 

In the scheme here considered, solar neutrino oscillations are generated by the mass- 
squared difference between v 2 and V\. It is clear from Eq. (|2.71 ) that the survival of solar f e 's 
mainly depends on the mixing angle $12, whereas the mixing angles $23 and i?24 determine 
the relative amount of transitions into sterile v s or active and v T . Let us remind the 
reader that and v T cannot be distinguished in solar neutrino experiments, since their 
matter potential and their interaction in the detectors are due only to neutral-current weak 
interactions and therefore they are equal. Thus, instead of and u T , one can consider the 
linear combinations 




cos d sin d 
— sin d cos d 




(2.8) 



with 



cosfi = ; 23 = , sin?9 = . 24 23 = . (2.9) 

1 r- 2 r- 2 n r 2 r 2 

L c 23 c 24 V c 23 c 24 



The mixing of v a and with V\ and v 2 is given by 



U a i = -S12V 1 - 44 , U a2 = C12 \jl - c| 3 c| 4 , U bl = U b2 = 0, (2.10) 
so that the solar neutrino oscillations occur only between the states 
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u e -»• u a with u a = c 23 c 24 v s + yl - c| 3 c| 4 ^ a , (2.11) 
with mixing angle $12, and 

44 = 1 - l^il 2 - \u a2 \ 2 = \u sl \ 2 + \U s2 \ 2 (2.12) 

gives the size of the projection of the sterile neutrino onto the state in which the solar v e 
oscillates. 

Therefore, solar neutrino oscillations depend only on -Qyi and the product C23C24. We 
distinguish the following limiting cases of Eq. (|2.11|) : 

• if C23C24 = then U s \ = U s2 = 0, U a \ = — s 12 , U a2 = Ci 2 , corresponding to the limit of 
pure two-generation v e — > v a transitions; 

• if C23C24 = 1 then U s \ = Si 2 , U s2 = C\ 2 and U a \ = U a2 = and we have the limit of 
pure two-generation v e — > v s transitions. 

Since the mixing of v e with v\ and v 2 is described by a two-generation equation with 
mixing angle $12, the mixing of v s with V\ and v 2 is given by the same formula multiplied 
by C23C24, and the mixing of v a with z/ x and v 2 is again described by the same equation 



times yl — C23C24, it is clear that, in the general case of simultaneous v e — > v s and v e — > v a 
oscillations, the corresponding probabilities are given by 

jjSun+Earth ^2 ^2 A p>Sun+Earth\ ic^ ^\ 

^ a Earth = (i - 44) (i - ^r rth ) , (2-i4) 

where P^^tf arth takes the standard two-neutrino oscillation form p4] for Am^ and ^12, 
found by numerically solving the evolution equation in the Sun and the Earth matter with 
the modified matter potential 

A = A cc + cl 3 c 2 24 A NC , (2.15) 

with Acc = 2\f2GpEN e and A NC — —V2GpEN n . Here N e and N n are, respectively, the 
number densities of electrons and neutrons in the medium, E is the neutrino energy and 
Gf is the Fermi constant. These expressions satisfy the relation of probability conservation 

pSun+Earth _i_ pSun+Earth i pSun+Earth i 

The analysis of the solar neutrino data in the four-neutrino mixing schemes is therefore 
equivalent to the two-neutrino analysis but taking into account that the parameter space 
is now three-dimensional (Am| 1; tan 2 $12, c^c^). Although originally this derivation was 



performed in the framework of the 2 + 2 schemes flT8| , |26|| , it is equally valid for the 3 + 1 
ones |37| . 

Let us finally comment on the range of variation of the mixing parameters. Choosing 
the convention that Am 2 ^ > 0, from Eqs. (|2.13|) and ( |2.14|) it is clear that for solar neutrino 



oscillations the full physical space is covered by choosing the mixing angles in the ranges: 

< ^12 < I , 0<1?23<^, 0<1?24<^. (2.16) 

Actually, since the dependence on $23 and $24 enters only through the combination c^c^, 
for what concerns the solar neutrino phenomenology one of these mixing angles could be 
fixed to zero. However, as we will see in the next section, they enter independently in the 
relevant probabilities for atmospheric neutrinos. 
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B. Atmospheric Neutrinos 



Atmospheric neutrino oscillations are generated by the mass-squared difference between 
is 3 and v±. From Eq. (|2.7|) one sees that the atmospheric v e 's are not affected by the four- 



neutrino oscillations in the approximation ^ 13 = i? 14 = and neglecting the effect of Am^ 
in the range of atmospheric neutrino energies. Conversely, the survival probability of at- 
mospheric z/^'s mainly depends on the mixing angle f? 34 . The mixing angles t?23 and "#24 
determine the relative amount of transitions into sterile v s or active v T . More explicitly 
Eq. (|2.7|) implies that the atmospheric neutrino oscillations, i.e. oscillations with the mass 



difference Am\ A and mixing angle $34, occur between the states 

up — > z/ 7 with up = s 23 C24^ s + c 2 3^ - s^s^v r and = s 24 ^ + c 24 ^r • ( 2 -17) 
We distinguish the following limiting cases of Eq. (|2 .17 ): 



• if C23 = 1 then U^i = = 0. The atmospheric = up state oscillates into a state 
^7 = C24^ T + S24^ s - This is the limit studied in Refs. [|^,|^]. We will denote this case 
as "restricted" . In particular: 

— in the case c 2 3 = c 2 4 = 1 we have U s3 = f/ s4 = 0, corresponding to the limit of 
pure two-generation — > v T transitions; 

— in the case c 2 3 = 1 and c 24 = we have U r3 = U T ± = 0, corresponding to the limit 
of pure two-generation — > v s transitions; 

• if C23 = there are no atmospheric neutrino oscillations as the projection of over 
the relevant states cancels out (U^ = = 0). 

Thus the mixing angle $23 determines the size of the projection of the over the "atmo- 
spheric" neutrino oscillating states, 

4 3 = Ki\ 2 + \u, 2 \ 2 = 1 - 1^ 3 | 2 - (2.18) 

As a consequence, one expects S23 to be small in order to explain the atmospheric neutrino 



deficit and, as we will see in Sec. IV, this is the case. Furthermore, the negative results 



from the CDHSW [34] and CCFR [35|] experiments on searches for ^-disappearance also 



constrain such a projection to be 

IM 2 + l^l 2 ^0.2 (2.19) 



"2:-! 



at 90% CL for Am[ SND > 0.4 eV 2 . 

Approximate expressions for the relevant survival probability including matter effects 
in the Earth can be found in Ref. []1|. For what concerns the work presented here, the 



relevant probabilities can be calculated by numerically integrating the evolution equation in 
the Earth with a modified matter potential, for example the diagonal piece of the potential 
takes the form 

A = (4 - sj 3 cl 4 )A NC , (2.20) 
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so that for pure atmospheric — > v T oscillations (s| 3 = = 0) A = 0, while for — > v s 
( s 23 = 0; s 24 = 1) -A = ^tvc*, as expected. It is this modification of the Earth matter 
potential that gives the atmospheric neutrino data the capability to discriminate between 
the active and sterile oscillation solution. In particular, higher sensitivity to this potential 
effect is achievable for the higher energy part of the atmospheric neutrino flux, which lead 
to the upward going muon data. The main effect of the presence of this potential is that for 
pure active oscillations the angular distribution of upgoing muons is expected to be steeper 
at larger arrival angles, while a flattening is expected from the matter effects for sterile 
oscillations |2"T| . 

From Eqs. ( |2.17| ) and Q2.20| ) we see that, unlike the case of solar neutrinos, the angles 
$23 and $24 enter independently in the atmospheric oscillations. Thus the analysis of the 
atmospheric neutrino data in the four-neutrino mixing schemes is equivalent to the two- 
neutrino analysis, but taking into account that the parameter space is now four- dimensional 

(Aml 3 ,$ u ,cl 3 ,4 A ). 

Concerning the range of variation of the mixing angles, the full parameter space relevant 
to atmospheric neutrino oscillation can be covered by choosing $23 and $24 in the ranges 
indicated in Eq. ( |2.16| ), but, in general, the mixing angle $34 must be allowed to take values 



in the interval 

7T 7T 



2 - J3A ~ 2 { J 



(in the convention Am\ z — 0)- I n the limiting case $23 = 0, there is an additional symmetry 
in the relevant probabilities so that the full parameter space can be spanned by < $34 < |, 
as expected since in this case we have the effective two-neutrino oscillation z/ M — > z/ 7 . 



III. SOLAR NEUTRINO ANALYSIS 

We first describe the results of the analysis of the solar neutrino data in terms of v e 
oscillations. We determine the allowed range of oscillation parameters using the total event 
rates of the Chlorine ||, Gallium |9HTT| as well as the recent SNO result on the CC event 
rate [T2\. For the Gallium experiments we have used the weighted average of the results 
from GALLEX+GNO and SAGE detectors. We have also included the Super-Kamiokande 
electron recoil energy spectrum measured separately during the day and night periods (cor- 
responding to the 1258-days data sample). This will be referred to in the following as the 
day-night spectra data, which contain 19 + 19 data bins. The analysis includes the latest 
standard solar model fluxes, the BP00 model [38], with updated distributions for neutrino 
production points and solar matter density. 

Details of the statistical analysis applied to the different observables can be found in 
Refs. 



In particular, we follow the approach of Ref. [25| and we do not include 



here the Super-Kamiokande total rate, since to a large extent the total rate is represented 
by the flux in each of the spectral energy bins. We define the x 2 function for the global 
analysis as: 



Xsol 



E 

i,i=l,41 



(R 



exp* 



G,ij 



(R 



th 
3 



Rf p ) 



(3.1) 
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where c/q ^ = a 2 R i j + Osp,ij- Here gr^ is the corresponding 41 x 41 error matrix containing 
the theoretical as well as the experimental statistical and systematic uncorrelated errors 
for the 41 rates while Gg p ,ij contains the assumed fully-correlated systematic errors for the 
38 x 38 submatrix corresponding to the Super-Kamiokande day-night spectrum data. We 
include here the energy independent systematic error which is usually quoted as part of the 
systematic error of the total rate. The error matrix Gr^ includes important correlations 
arising from the theoretical errors of the solar neutrino fluxes, or equivalently of the solar 
model parameters. Notice that this statistical treatment is not the same employed in our 
pre-SNO analysis |JT] what makes the comparison not totally straight forward (see Ref. 
for a details on the comparisons). 

As discussed in the previous section, the analysis of the solar neutrino data in these four- 
neutrino mixing schemes is equivalent to the two-neutrino analysis but taking into account 
that the parameter space is now three-dimensional (Am| 1; tan 2 $ 12 , c^c 2 ,^), so that we have 
a total of 38 degrees of freedom (d.o.f.). 

We first present the results of the allowed regions in the three-parameter space for the 
global combination of solar observables. Note that, since the parameter space is three- 
dimensional, the allowed regions for a given CL are defined as the set of points satisfying 
the condition 

xU^™ 2 2 , tf 12 , 44) - X 2 ol , mi n < A X 2 (CL, 3 d.o.f.) (3.2) 

where, for instance, A^ 2 (CL, 3 d.o.f.) = 6.25, 7.81, and 11.34 for CL = 90%, 95% and 99% 
respectively, and Xsol min * s the global minimum in the three-dimensional space. 

In Fig. |2| we plot the sections of such a volume in the plane (Ara^, tan 2 ($12)) for different 
values of the active-sterile admixture \U s i\ 2 + \U s2 \ 2 = c^c 2 ,^. The global minimum used in 
the construction of the regions lies in the LMA region and for pure v e -active oscillations 

\U sl \ 2 + \U s2 \ 2 = cj 3 c 2 24 = 

Am 2 21 = 4.1 x 1(T 5 eV 2 
tan 2 #21 = 0.41 

Xsol.min = 35.3 (3.3) 

which for 38 d.o.f. (3 rates + 38 spectrum points — 3 parameters) corresponds to a goodness 



of the fit (GOF) of 59%. Notice that we have used the Chooz reactor bound to cutoff 
the allowed solutions for mass values above ~ (7-8) xlO~ 4 eV 2 . The first and last panels 
shown in Fig. correspond to the limiting cases of pure active and pure sterile oscillations 
presented in Fig.l of Ref. p5| . 

As seen in Fig. ^| the SMA region is still a valid solution for any value of c| 3 c| 4 at 99% 
CL. This is expected as in the two-neutrino oscillation picture this solution holds for both 
pure v e -active and pure f e -sterile oscillations. Notice, however, that the statistical analysis 
is different: in the two-neutrino picture the pure z/ e -active and v e -sterile cases are analysed 
separately, whereas in the four-neutrino picture they are taken into account simultaneously 
in a consistent scheme. On the other hand, both the LMA and LOW solutions disappear 
for larger values of the active-sterile admixture. 

Let us comment here that, unlike the case of atmospheric neutrinos, where the dis- 
crimination between active and sterile oscillations arises from the difference in the matter 
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potentials, for solar neutrinos the main source of difference is due to the lack of neutral- 
current contribution to the water Cerenkov experiments for the sterile case. Unlike active 
neutrinos, which lead to events in the Super-Kamiokande detector by interacting via NC 
with the electrons, sterile neutrinos do not contribute to the Super-Kamiokande event rates. 
Therefore a larger survival probability for 8 B neutrinos is needed to accommodate the mea- 
sured rate. As a consequence, a larger contribution from 8 B neutrinos to the Chlorine and 
Gallium experiments is expected, so that the small measured rate in Chlorine can only be 
accommodated if no 7 Be neutrinos are present in the flux. This is only possible in the SMA 
solution region, since in the LMA and LOW regions the suppression of 7 Be neutrinos is not 



sufficient. Comparing with our pre-SNO results in Ref. |31[] we see that, after the inclusion 
of the SNO result, which observes a CC event rate for 8 B neutrinos lower than the total 
event rate of Super-Kamiokande, the SMA solution either active or sterile is now less likely 
since it is difficult to fit simultaneously the Super-Kamiokande flat spectrum and their total 
event rate and the SNO CC measurement. 

Notice also that the SMA region for oscillations into sterile neutrinos is slightly shifted 
downwards with respect to the active case. This is due to the small modification of the neu- 
trino survival probability induced by the different matter potentials. The matter potential 
for sterile neutrinos is smaller than for active neutrinos due to the negative NC contribution, 
proportional to the neutron abundance. For this reason the resonant condition for sterile 
neutrinos is achieved at lower Am 2 . 

In Fig. [3] we plot the difference Ax 2 ol between the local minimum of \lo\ f° r eacn solution 
and the global minimum in the three-dimensional space, as a function of the active-sterile 
admixture \U s i\ 2 + \U S 2\ 2 = 44- From the figure we find that: 

• solar neutrino data favour pure u e — > v a oscillations but sizeable active-sterile admix- 
tures are still allowed; 

• the dependence of Ax 2 ol on the active-sterile admixture is very gentle in the SMA 
region while it is much more pronounced for the LMA solution; 

• the three-dimensional regions are acceptable at 90% (99%) CL for the following values 
of c 23 c 2 4 

44 < - (I-) for SMA, 

44 < 0.45 (0.72) for LMA, (3.4) 
c 2 3 c 2 4 < 0.40 (0.99) for LOW-QVO; 

• at 99% CL all the three-dimensional regions are allowed for maximal active-sterile 
mixing c 2 3 c 2 4 = 0.5. 

There is a subtlety on the meaning of these limiting values of the mixing c^c 2 ^. The values 
quoted in Eq. ( |3.4|) correspond to the mixings at which the corresponding three-dimensional 
region disappears at a given CL. Strictly speaking this is not the same as the allowed range of 
44 a ^ a given CL when the parameters (tan 2 $12, Am^) are left free to vary, either in the 
full parameter space or within a given solution. Such limits would be obtained by imposing 
the condition Ax s 2 ol = Xsoi.mm i 2 (44) ~ Xsoi.mm < Ax 2 (CL, 1 d.o.f.) where X s 2 i,mm i 2 (44) 
is minimized in the parameter space 1-2 (or in a given region of this parameter space) and 
Xsoimm i s minimized in the same region of the space 1-2 and in the mixing c^c^. 
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IV. ATMOSPHERIC NEUTRINO ANALYSIS 



In our statistical analysis of the atmospheric neutrino events we use all the samples of 
Super-Kamiokande data: e-like and ^x-like samples of sub- and multi-GeV || data, each 
given as a 5-bin zenith-angle distribution^, and upgoing muon data including the stopping 
(5 bins in zenith angle) and through-going (10 angular bins) muon fluxes. We have also 
included the latest MACRO [[32] upgoing muon samples, with 10 angular bins, which is also 
sensitive to the active-sterile admixture. So we have a total of 45 independent inputs. 

For details on the statistical analysis applied to the different observables, we refer to 



Refs. [[Hynj. As discussed in the previous section the analysis of the atmospheric neutrino 
data in these four-neutrino mixing schemes is equivalent to the two-neutrino analysis, but 
taking into account that the parameter space is now four- dimensional (Am| 3 , $34, c| 3 , c 2 ,^). 
Alternatively we present the results in the equivalent parameter space (Ara^,^, |^i| 2 + 
|^ 2 | 2 ,|^i| 2 + |C/ s2 | 2 ). 

We first present the results of the allowed regions in the four-parameter space for the 
global combination of atmospheric observables. Notice that since the parameter space is 
four- dimensional the allowed regions for a given CL are defined as the set of points satisfying 
the condition 

XatmlAm^, tf 34 , cj 3 , 4) - Xatm.min < ^(CL, 4 d.O.f.) (4.1) 

where Ax 2 (CL, 4 d.o.f.) = 7.78, 9.49, and 13.3 for CL = 90%, 95% and 99%, respectively, 
and Xatm min i s the global minimum in the four- dimensional space. 

In Fig. f| we plot the sections of such a volume in the plane (Am 43 , sin 2 (-i9 34 )) for different 
values of the mixings $23 and $24, which we parametrize by values of the projections | C/^i | 2 — I— 
I [/"^ 1 2 and \U s i\ 2 + \U S 2\ 2 - As discussed in Sec. [TT] for arbitrary values of c 2 3 and c 2 4 the full 
parameter space is covered for ?9 34 in the range given in Eq. ( J2.21Q . In Fig. |] we display 
this full parameter space by showing the region in (Am 43 , sgn($ 34 ) sin 2 ($ 34 )) with "positive" 
sin 2 ($ 34 ) for < $34 < | and "negative" sin 2 ($34) for — | < -# 34 < 0. As seen from the figure 
the parameter space is composed of two separated regions, each of them around maximal 
mixing $34 = ±f • 

The global minimum used in the construction of the regions lies almost in the pure 
atmospheric v^-Vt oscillations. At the best=fit point 

i^ir+if/ s2 r = c 2 3 c 2 4 = o.97 

|^i| 2 + |^2| 2 = 4 = 0.01 

Am 2 3 = 2.4 x 10~ 3 eV 2 
^34 = 39° 

Xatm.min = 29.0 (4-2) 

which for 45 — 4 = 41 d.o.f. corresponds to a GOF of 92%. A careful reader may find this 
GOF surprisingly high. Let us comment here that there has been a reduction of the x 2 



*Note that for convenience and maximal statistical significance we prefer to bin the Super- 
Kamiokande contained event data in 5, instead of 10 bins. 
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with respect to previous analyses |]19| , p9| due to a better agreement of the 79.5 SK electron 



distribution with the non-oscillation expectation (the same effect is found in Ref. [p9||). 
Since the scenarios discussed here do not affect the atmospheric v e predictions, this implies 
an overall improvement on the GOF of any of these four-neutrino oscillation schemes. 

From Fig. [| we see that the region becomes considerably smaller for increasing values 
of the mixing angle $23, which determines the size of the projection of the over the 
"atmospheric" neutrino oscillating states, and for increasing values of the mixing angle $24, 
which determines the active-sterile admixture in which the "almost-z/ M " oscillates. Therefore 
from the analysis of the atmospheric neutrino data we obtain an upper bound on both 
mixings, which, in particular, implies a lower bound on the combination c| 3 c| 4 = 
|f/ S 2| 2 limited from above by the solar neutrino data. To quantify these bounds we display 
in Fig. |]a the allowed region in the parameter space (s| 3 = |t/ M i| 2 + {U^ \ 2 , <^3 c m = Wsi\ 2 + 
\U S 2\ 2 )- Following the discussion below Eq. (|3.4j) we obtain the allowed two-dimensional 
region by imposing the condition 

A Xatm = Xatm.min 34^23, ^4) " xLn < A X 2 (CL, 2 d.O.f.) , (4.3) 

where x 2 tm mm 34(^231 ^24) is minimized with respect to Am 2 3 and $34, and X 2 tmmin ^ s the 
global minimum in the four-dimensional parameter space. The diagonal cut in Fig. |6|a gives 
the unitarity bound 

im 2 + i^r + i^ii 2 + pm 2 = i^ t3 i 2 + i^ 4 i 2 < 1 , (4.4) 

where we have used that in Eq. ( [2.7|) U e ^ = [7 e4 = 0. As seen from Fig. ||a, the analysis 
of the atmospheric data imposes a severe bound on the projection. For instance, the 
two-dimensional region extends only to 

4 = l^ir + l^2| 2 ^0.12(0.16) (4.5) 

at 90% (99%) CL, which, as previously mentioned (see Eq. (|2.19|) ), is perfectly consis- 
tent with the bounds from the disappearance accelerator experiments CDHDW @] and 
CCFR |H given in Eq. (p39j ). 



In order to determine the impact of the small but possible deviation from zero of s\\ we 
also study the "restricted" case of Refs. p8| , p9| in which we impose s\\ = \U^i\ 2 + \U^2\ 2 = 0. 
In this case the parameter space for atmospheric neutrino oscillation is three-dimensional 
and the allowed regions for a given CL are defined as the set of points satisfying the condition 

Xatm( A ™L ^34, 4^) - xltm,min < A X 2 (CL, 3 d.O.f.) . (4.6) 

In Fig. [| we plot the sections of this volume in the plane (Am 2 3 , sin 2 ($34)) for different values 
of c| 4 . As discussed in Sec. [Tl] for the "restricted" case the full parameter space is covered 
with < $34 < |. This is graphically displayed in Fig. where the symmetry "#34 — > — -$34 
is evident. 

The constraints on the active-sterile admixture in the two cases are displayed in 
Fig. ||b, where we plot the dependence of Axa tm 011 the active-sterile admixture c 2 3 C2 4 = 
|^si| 2 + |^s2| 2 , both in the general case (i.e. when xttm * s minimized with respect to all other 
parameters in the problem) and in the "restricted" case (for which the projection has 
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been fixed \U^i\ 2 + {U^l 2 = 0). The horizontal lines are the values of Ax 2 (CL, 1 d.o.f.) for 
CL = 90% and 99%. From this figure we obtain the 90% (99%) CL lower bounds on d^c^ 
from the analysis of the atmospheric neutrino data: 

c\Ai = \Usi? + \U s2 \ 2 > 0.64 (0.52) (4.7) 
c 2 4 = \U sl \ 2 + \U s2 \ 2 > 0.83 (0.74) for the restricted case. (4.8) 

As seen in the figure Ax 2 tm nas a monotonically growing behaviour as |f/ s i| 2 + |L r S 2| 2 decreases 
for |l/si| 2 + \U S 2\ 2 ^ 0.4-0.6, while the dependence flattens below that value with the presence 
of a shallow secondary minimum around \U s i\ 2 + \U S 2\ 2 ~ 0.25-0.3. We have traced the 
presence of this secondary minimum to the combined effect of the behaviour of x 2 f° r the 
different sets of upgoing muon data. For all sets of upgoing muon data the best-fit is 
very close to the global minimum in Eq. ([4.2]) . MACRO data at the moment lead to a 
monotonically growing % 2 , as the sterile admixture grows while Super-Kamiokande upgoing 
data on both stopping and through-going muons presents a flattening of their \ 2 f°r |£^i| 2 + 
\U S 2\ 2 < 0.5-0.6. This loss of sensitivity of Super-Kamiokande at larger sterile admixtures 
arises from the fact that their latest upgoing muon data presents a flattening of the muon 
suppression in the angular bins for arrival directions around —0.7 < cos# < —0.4, which is 
compatible with large sterile admixtures. The same effect has been found in Ref. ||29| . When 
combining the x 2 f°r MACRO and Super-Kamiokande data samples, and taking into account 
the strong correlation between their respective theoretical errors, we find the appearance of 
this secondary shallow minimum. 

From the figure we see that oscillations into a pure sterile state are excluded at 99.97% 
CL as compared to pure active oscillations. One must notice, however, that taken as an 
independent scenario, — > v a oscillations give a Xmin = 46 (for 45 — 2 d.o.f.), which implies 
a GOF of 35% and therefore cannot be ruled out on the basis of its absolute quality. 

In summary, we see that the analysis of the atmospheric data implies that: 

• one of the two- neutrino states in the 3-4 oscillating pair must be a close-to-pure 

• this close-to-pure oscillates into a state which is preferably composed by p T , although 
a 48% admixture of sterile in this state is still allowed at 99% CL; 

• when imposing the "restricted" condition that one of the states is a pure u^, the bound 
on the sterile admixture is tightened approximately by a factor 2 (from 48% to 26% 
at 99% CL). 



V. COMBINED ANALYSIS: CONCLUSIONS 

From the previous sections we have learned that the analysis of the solar data favours 
the scenario in which the solar oscillations in the plane 1-2 are v e oscillations into an active 
neutrino, and from that analysis we found an upper limit on the projection of the v s on the 
1-2 states. On the other hand, the atmospheric neutrino analysis prefers the oscillations of 
the 3-4 states to occur between a close-to-pure and an active (v T ) neutrino, thus giving 
an upper bound on the projection of the v s over the 3-4 states, or equivalently a lower bound 
on its projection over the 1-2 states. The open question is then what the best scenario is 
for the active-sterile admixture once these two bounds are put together. 
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To address this question in full generality we have studied the behaviour of the global 
Xatm+soi function defined as 

Xatm+ S ol( Am 21> #12, A ™L # 34 , # 23 , #24) = 

Xsol( A "4^23, £23^4) + Xatm( Am 4 3l ^34, 4 3 , c\ A ) (5.1) 

with the mixing angles $23 and $24 determining the active-sterile admixture for the following 
scenarios: 

• ATM+SOL L ma(sma) : in this case we minimize X 2 tm+sol with respect to Am^, i?i 2 , 
Am^ an d #34 with Am^ and $12 constrained to be in the LMA (SMA) region. The 
Up projection S23 — l^i| 2 + I ^2 1 2 is f ree to vary in the range allowed by the analysis. 

• ATM r +SOL LMA (sma): in this case we fix S23 = |^i| 2 + |^2| 2 = 0, so \U sl \ 2 + \U s2 \ 2 = 
C24. We minimize X 2 tm+so i with respect to Am^, $12, Am^ and $34, with Am^ and 
f?i2 constrained to be in the LMA (SMA) region. 

In Table | we give the values Xminatm+soi °f the best-fit point and its associated GOF in 
the six-dimensional (five-dimensional for the "restricted" case) space in each of these cases 
together with the best-fit value of the mixings |J7 M i| 2 +|£^t2| 2 = and |f/ s i| 2 +|t/ S 2| 2 = C23 c 24- 
In Fig. |7] we plot the shift in \ 2 from the corresponding global minimum in each of the 
scenarios as the mixing \U s i\ 2 + \ U S 2\ 2 = c\ z c\ A varies. From the table and the figure we find 
the following behaviours: 

• There are two favourite configurations, which we will denote as near-pure-sterile 
solar neutrino oscillations plus near-pure-active atmospheric neutrino oscillations 
(NPSS+NPAA), and close-to-active solar neutrino oscillations plus close-to-sterile at- 
mospheric neutrino oscillations (CAS+CSA), respectively. NPSS+NPAA oscillations 
are characterized by the minima 

\U sl \ 2 + \U s2 \ 2 ~ 0.91-0.97, (5.2) 

where the exact position of the minimum depends on the allowed admixture |£/ M i| 2 + 
If/^l 2 - CAS+CSA oscillations are characterized by the minima 

\U s i\ 2 + \U s2 \ 2 ~ 0.18-0.2 . (5.3) 

In both cases the minima are not very deep in x 2 (this is particularly the case for 
CAS+CSA). 

• In none of the cases is the maximal active-sterile (MAS) admixture \U s i\ 2 + \U S 2\ 2 = 
c 23 c 24 = 0-5 favoured. 

The scenarios with the solar solution within the LMA region (Fig. |7|b) prefer the CAS+CSA 
scenario, although the dependence of X 2 tm+so i with the active-sterile admixture is not very 
strong and it presents a secondary minimum near the complementary situation of close- 
to-sterile solar plus close-to-active atmospheric (CSS+CAA) scenario, which is acceptable 
at 90% CL. They give the best fit to the combined analysis with a GOF: 67% for the 
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"unrestricted" case. The "restricted" case gives a worse fit but still with an acceptable GOF 
of 59%. In these scenarios MAS admixture is ruled out at 99% CL for the "restricted" case 
and it is acceptable at about 95% CL for the unrestricted one. This behaviour arises from 
the fact that in these scenarios increasing the sterile admixture both in the atmospheric pair 
and in the solar case is strongly disfavoured. It is precisely in this case that one would expect 
the combined solution to lie near MAS admixture. However, the results show that it is the 
solar dependence (which dominates in the combination) that leads to the slightly better 
description for the CAS+CSA scenario with a secondary minimum at the complementary 
CSS+CAA scenario. 

On the other hand, scenarios with the solar solution within the SMA region (Fig. 0b) 
prefer the NPSS+NPAA scenario. They give the worse fit to the combined analysis for the 
"unrestricted" case with a GOF of 62% improving to 65% for the "restricted" case . In these 
scenarios the dependence of Xatm+soi with the active-sterile admixture is more pronounced 
and it is dominated by the dependence of the atmospheric neutrino analysis, although there 
is a secondary minimum at the CAS+CSA configuration, which for the '"unrestricted" case 
is acceptable at 99% CL. In all these scenarios MAS admixture is ruled out at more than 99% 
CL. Qualitatively this behaviour arises from the fact that, in these scenarios, increasing the 
sterile admixture in the atmospheric pair is strongly disfavoured while it is still acceptable 
for the solar pair since the SMA solution is valid for pure v e — > v s solar oscillations. 
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FIG. 1. 2 + 2 mass schemes in four- neutrino mixing favoured from the combination of solar, 
atmospheric and LSND results 
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ATM+SOL L ma 


0.065 


0.21 


73.8 


67% 


ATM r +SOL L ma 
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0.22 


77.4 


59% 


ATM+SOLsma 


0.030 


0.91 


75.4 


62% 


ATM r +SOL S ma 


0. 
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75.5 
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TABLE I. Xmin atm+sol °f the best-fit point and its associated GOF in the six-dimensional 
(five-dimensional for the "restricted" case) space in each of the cases discussed in the text together 
with the best fit value of the mixings s| 3 = |t^ti| 2 + I 1 2 and c^c 2 ,^ = \U s i\ 2 + \U S 2\ 2 - 
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FIG. 2. Results of the global analysis of solar neutrino data for the allowed regions in Am|i and 
tan 2 $12 for the four-neutrino oscillations. The different panels represent sections at a given value 
of the active-sterile admixture \U S \\ 2 + \ U S 2\ 2 = c^c^ of the three-dimensional allowed regions at 
90%, 95% and 99% CL. The best-fit point in the three-parameter space is plotted as a star. 
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FIG. 3. Difference A% 2 between the local minimum of x 2 f° r each solution of the global analysis 
of solar neutrino data and the global minimum in the plane (tan 2 $12, Ara^) as a function of the 
active-sterile admixture |f7 s i| 2 + \U S 2\ 2 = c\nC^ 4 . 
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FIG. 4. Results of the analysis of atmospheric neutrino data for the allowed regions in Am| 3 
and $34 for the four- neutrino oscillations. The different panels represent sections at given values of 
the projection If/^iP + lt/^l 2 = «23 an d the- active-sterile admixture |t/ s i| 2 + |?7s2| 2 = c 23 c 24 °f the 
four-dimensional allowed regions at 90%, 95% and 99% CL. The best-fit point in the four-parameter 
space is plotted as a star. The last panel corresponds to the case in which x 2 has also been 
minimized with respect to s 2 .^ and c^c 2 .^. 
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FIG. 5. Results of the analysis of atmospheric neutrino data for the allowed regions in 
Am 2 ^ an d "$34 for the four-neutrino oscillations in the "restricted" case of projection 
Wni 1 2 + 1^2 1 2 = = 0. The different panels represent sections at given values of the ac- 
tive-sterile admixture \U S \\ 2 + \ U S 2\ 2 = c| 4 of the three-dimensional allowed regions at 
and 99% CL. The best-fit point in the three-parameter space is plotted as a star. 
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FIG. 6. (a) 90%, 95% and 99% CL allowed regions for the projection |C/ M i| 2 + \U^ 2 \ 2 = s 2 3 
and the active-sterile admixture |t4i| 2 + |t/ S 2| 2 = c 23 c 24 from the analysis of the atmospheric data 
(see text for details), (b) Ax 2 as a function of the active-sterile admixture \U s i\ 2 + \U S 2\ 2 = C23C24. 
In the lower curve x 2 nas been minimized with respect to $34, Am| 3 and the Vu projection 
|C^ii| 2 +|t^i2| 2 = S23- The upper curve corresponds to the "restricted" case | t/^i | 2 + 1 f/^2 1 2 = s 2 3 = 0, 
for which x 2 has been minimized only with respect to $34 and Am 2 3 . 



21 




FIG. 7. A X 2 as a function of the active— sterile admixture |t/ s i| 2 + |C/ S 2| 2 — c 23 c 24' 

In 

the ATM curves x 2 have been minimized with respect to $34, Am^ and the Vu projection 
|£^i| 2 + I Ua2 1 2 = s 2 3, as well as with respect to the 1-2 parameters $12 and Amfi- The ATMr 
curves correspond to the "restricted" case |?7«i| 2 + |fu2P = s 23 = 0> so X 2 has been minimized only 
with respect to #34, Am 2 3 and the 1-2 parameters. In the upper (lower) panel the 1-2 parameters 
$12 and Am^i have been constrained to lie in the LMA (SMA) solution region. 
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